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Abstract. The evolution pattern of exceptional points is studied in a non-integrable limit 
of the complex-extended 3-level Richardson-Gaudin model. The appearance of a pseudo- 
diabolic point from the fusion of two exceptional points is demonstrated in the anti-hermitian 
limit of the model and studied in some details. 



1 Introduction 

Small quantum systems, whose properties are profoundly affected by environment, 
i.e., continuum of scattering and decay channels, are intensely studied in vari- 
ous fields of physics (nuclear physics, atomic and molecular physics, nanoscience, 
quantum optics, etc.). These different open quantum systems (OQS), in spite of 
their specific features, have generic properties, which are common to all weakly 
bound/unbound systems close to the threshold. An essential part of the motion of 
short-lived nucleonic matter is in classically forbidden regions, and their properties 
are impacted by both the continuum and many -body correlations [1,2]. 

Resonances are commonly found in various quantum systems, independently 
of their building blocks. Resonances are genuine intrinsic properties of quantum 
systems describing preferential decays of unbound states. The effect of resonances 
and the non-resonant scattering states can be considered in the OQS extension of the 
shell model (SM), the so-called continuum shell model (CSM) [2]. Two realizations 
of the CSM have been studied recently: the real-energy CSM [2-A] and the complex- 
energy CSM [6,7] based on the Berggren ensemble [5], the so-called Gamow Shell 
Model (GSM). 

For hermitian Hamilton operators, both real-energy CSM [2-A] and complex- 
energy CSM (GSM) [6,7] lead naturally to the non-hermitian (complex-symmetric) 
eigenvalue problem. As a result, OQSs exhibit several unintuitive properties, which 
make them qualitatively different from closed quantum systems (CQS). Among 
those salient features are: the segregation of time scales in the continuum [11, 12] 
(see also Ref. [2] for a recent review), the alignment of near-threshold states with 
decay channels [1,2, 13], the instability of SM eigenstates at the channel thresh- 
old [14, 15], or the resonance crossings [2, 16]. In this lecture, we shall concen- 
trate on the latter phenomenon which will be illustrated on the example of a non- 
integrable pairing model. 
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Most studies of degeneracies associated with avoided crossings in quantal spec- 
tra, focused on the topological structure of the Hilbert space and the geometric 
phases [17, 18]. Among these degeneracies, one finds a diabolic point (DP) in her- 
mitian Hamiltonians [17, 19], and an exceptional point (EP) [20-22], which appears 
in the complex g-plane of a generic Hamiltonian H(g) = H o + gH\, where both 
Ho and Hi are hermitian and [Ha, Hi] ^ 0. Below, we shall introduce a prototyp- 
ical OQS, the 3-level Richardson-Gaudin (RG) model, to discuss the appearance of 
resonance crossings and their evolution with a parameter of a non-integrable pertur- 
bation. In particular, we shall show the appearance of a new kind of degeneracy, a 
pseudo-diabolic point (pseudo-DP) in the anti- hermitian limit of this model. 



2 The 3-level pairing model 

RG models [23] are based on the ST/ (2) algebra with elements K~f, , and Kf, 
fulfilling the commutation relations: [K+.K^] = 5 W K? , [K?,K±] = ±5 w Kj t , 
where indices /, /' refer to a particular copy from a set of L, SU(2) algebras. Each 
SU (2) algebra possesses one quantum degree of freedom. In the following, we shall 
use the pair representation of SU(2) algebra leading to pairing Hamiltonians. The 
elementary operators in this representation are the number operators Nj and the pair 
operators Aj, Aj, defined as: 

N i = 2 a Uajm ; At = a] m a] m = (At)t (1) 

m m 

where j is the total angular momentum and m is the z-projection. The state jm is 
the time reversal of jm. The relation between the operators of the pair algebra and 
the generators of the SU(2) algebra is: 

K? = ±N l -\f2 l ; K+ = (K^ = \A\ (2) 

where fii is the particle degeneracy of level I. With this correspondence, one can 
introduce an integrable 3-level pairing Hamiltonian as: 

ff/(9)=^^ i +S^A,U j (3) 

i ij 

Below, we shall consider the non-integrable version of the 3-level pairing model 
[24]: 

H NI {g) = H I {g)+g'Y,rf , g' = 75 (4) 



where 7 < is a C-number. 
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Figure 1. Localization of EPs (the red cross) and pseudo-DPs (the red cross inside a blue 
circle) in the complex (?-plane for the non-integrable 3-level pairing Hamiltonian with g' / g = 
— 1/2 (see Eq. (4)). For chosen energies of levels, one finds a trivial degeneracy at g — (the 
non-interacting limit). For more details, see the description in the text. 



3 Level degeneracies 

The position of all possible degeneracies in the complex g-plane are indicated by 
the roots of the coupled equations: 

det[H(g)-EI}=0; ^det [H (g) — EI] = (5) 

By eliminating E from these two equations, we are left with the discriminant D(g), 
a polynomial in g of degree M = n(n — 1), where n is the number of eigenstates. 
The discriminant can be written as [21]: 

D(g) = J] [E m (g)-E m ,(g)} 2 (6) 

m<m' 

where E m (g), E m >(g) denote the complex eigenvalues of H(g). The eigenvalue 
degeneracies E m (g) = E m >(g) at g — g a {a — 1, . . . , M), can be found numer- 
ically by looking for sharp minima of the functional D(g). For the non-integrable 
Hamiltonian, the degenerate eigenvalues are either the single-root (EP) or double- 
root solutions such as DP or pseudo-DP. 
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Figure 2. The formation and decay of a pseudo-DP at g' = —g/2 (the blue star) in the 
anti-hermitian limit of the non-integrable 3-level pairing Hamiltonian (4). 

4 Properties of a pseudo-diabolic point in the anti-hermitian 
limit of a non-integrable 3-level pairing model 

Let us solve the non-integrable 3-level pairing model for the case of 2 pairs of 
fermions in a valence space of degeneracy Q\ = 2, i?2 = 6, fij, = 2. The ener- 
gies of three levels are e% = 0, £2 — 1, and £3 = 2. In this model space, there are 4 
many-body states. 

Fig. 1 shows the global pattern of level crossings in the complex g-plane for 
g'/g = —1/2. In the lower half-plane of g, all eigenvalues are either discrete states 
at the real g-axis or decaying resonances. Complex conjugate degeneracies situated 
in the upper half-plane (Im(g) > 0) correspond to capturing resonances. 

Fig. 2 shows a typical pattern associated with the formation/breakup of the 
pseudo-DP at g = — i/(4v / 2). With an increasing value of g'/g, the two EPs in 
the complex g-plane approach each other, coalesce at g'/g = —1/2 and, subse- 
quently, move along the axis IZe(g) = for g'/g > —1/2. The coalescence of two 
EPs at g'/g = — 1/2, leads to the formation of a double-root singularity, for which 
the geometric phase associated with a cyclic evolution changes twice faster than for 
the EP. On the other hand and in contrast to the ordinary DP, both eigenvalues and 
eigenvectors merge, like for an ordinary EP. 

The upper part of Fig. 3 shows real parts of phase trajectories around the pseudo- 
DP for all four eigenvectors in the model space. At the pseudo-DP, real part of 
their eigenenergies are identical. These trajectories in g-plane are described by: 
g = — (l/(4\/2)i + 0.01e 1¥> , where tp defines the way of encircling the degeneracy. 
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Figure 3. The real parts of phase trajectories for four eigenvectors involved in the pseudo-DP 
at 9 /9 = — 1/2 ( me upper part). This degeneracy of two eigenvalues results from a coales- 
cence of two EPs. In the lower part, the real parts of phase trajectories for two eigenvectors 
which form an EP at g'/g = —49/100 (see Fig. 2) are shown. The phase of each vector is 
defined with respect of its reference value at ip = 0. 



The phase of each vector is defined with respect of its value at (p = 0. The eigenval- 
ues for vectors '2' and '3' are degenerate at the pseudo-DP and change their phases 
by 7r at each encircling. Simultaneously, phases of vectors for non-degenerate eigen- 
values '1' and '4' remain approximately constant. For a comparison, the lower part 
of Fig. 3 exhibits real parts of phase trajectories of all four vectors around the EP 
g = — 0.207687i at g'/g = —49/100 (see Fig. 2). These trajectories are specified 
by: g = — 0.207687i + 0.01e lv . Eigenvalues for vectors '2' and'3' are degenerate at 
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Figure 4. Real and imaginary parts of eigenvalues at around the pseudo-DP g — — i/ (4\/2) 
are plotted along the cut {JZe(g) 1 — l/(4\/2)) in g-plane. For more details, see the description 
in the text. 



the EP. After first encircling of the EP they become exchanged. After a second encir- 
cling, the phases of vectors '2' and '3' change by it, and after two more encircling 
of the EP all phases return to their initial values. 

The behavior of eigenenergies in the neighborhood of a pseudo-DP in complex 
g-plane are shown along the cut (IZe(g), — 1/(4\/2)) (see Fig. 4). The fusion of 
two EPs and the formation of a single pseudo-DP is seen as a sharp crossing of four 
energies and a coalescence of two widths (the eigenvectors '2' and '3'). Even though 
phases of eigenvectors '1' and '4' are almost constant when encircling the pseudo- 
DP (see Fig. 3), they are essential for the formation of a pseudo-DP degeneracy. It 




can be shown rigorously, that the pseudo-DP does not result from the Jordan block 
of 2-level system. 



5 Mixing of wave functions at the pseudo-diabolic point: the 
pairing correlation energy 

Salient features of eigenvectors around a pseudo-DP of a 3-level pairing Hamilto- 
nian (4) with g'/g = — 1/2 can be studied analytically. Let g — — i(4v / 2) _1 + 8, 
where 8 is a complex number. The eigenenergies of a Hamiltonian matrix: 



E 2 
E 3 

w 



35.9338 8 + 0(8)' 
0{8f 

^2 



/4- 3. 79878 i 

4- %/2i + 8<5 
4- V2i + 0(S) 
0.263243 i - 7.93378 8 + 0{8) 2 J 



\4 



(7) 



depend parametrically on the parameter 8. Notice an asymmetric dependence on 8 
of eigenvalues '2' and '3' which form a pseudo-DP. 

The eigenvectors, which are normalized according to the dual metric of RHS 
[5,8-10]: 



< ili\Uj >= Si 



(8) 



can be written as: 
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0.616894 + 0.517406 i+ (7.05076- 1.17993 i) <5 - (27.8062 + 130.27 i) S 2 + O(S) 3 ' 1 
0.731723 - 5.83203 i 5 - 127.824 S 2 + 0{S) 3 
0.488757 - 3.21838 i 5 - 67.7276 S 2 + 0(8) 3 

.0.616894- 0.517406 i- (7.05076+ 1.17993 i) S - (27.8062 - 130.27 i) S 2 + 0(8f 



u 2 



(l-i) 
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4V2\/2 V* 8V2V2 



l-i _ 

"42T74 



7? + ^±=V^+0(^) 5 

4V2\/2 



M4 



^371 (1 + i) Vtf + H \/^# (1 - i) J* + 0(<5) ! 

-J4L. 1 + 9 y^+8- (9V2-S) i ^ + 3 
4V2V2 v* 8V2V2 

0.345603 - 0.412057 i- (2.62262- 1.08314 i) <S + (14.2829 + 10.4398 i) (5 2 + 0(5) E 
-0.582736 + 0.412499 i <5 - 6.21696 S 2 + 0(8) 3 
-0.87242 + 3.06004 i 5 + 26.4056 S 2 + O(S) 3 
.0.345603+ 0.412057 i+ (2.62262+ 1.08314 i) S + (14.2829 - 10.4398 i) S 2 + 0(5f 



In the limit 5 — > 0, the eigenvectors M2 and {43 with (identical) eigenvalues 
i? 2 , -E3 (see Eq. (7)) have divergent components. This feature of eigenvectors at a 
pseudo-DP makes it similar to an ordinary EP and leads to a singular behavior of 
various quantities. 

Let us consider for example the pairing operator: g J2ij A\Aj, in the neighbor- 
hood of a pseudo-DP for g'/g = — 1/2: 

(a 2 i + b 2 6 + O(S) 2 a 5 /VS + O(S) 1 / 2 a 6 /VS + O(S) 1 / 2 a 4 i + b 4 S + 0(S) 2 \ 
a 5 /VS + 0(5y/ 2 ai /5 + 0(6) ai i/S + 0(6) a 7 /VS + O(S) 1 / 2 
ae/VS + OiS) 1 / 2 aii/5 + 0{5) - ai /5 + 0{8) a s /VS + OiS) 1 ' 2 
Va 4 i + ^ 5 + O(S) 2 a 7 /V5 + 0(5)^ 2 a 8 /VS + 0(5)^ 2 a 3 i + b 3 S + 0(8) 2 J 

Coefficients Oj in this expression are: 

ai = 4:, »2 = -7.43796, 6 2 = -3.64018, 
16 

a 3 = 0.455281, 63 = 0.765176, a 4 = 0.603023, 64 = 1.8219, 

a 5 = 0.475579 (1 - i), 65 = 0.449184 (1 + i), a 6 = 0.475579 (1 + i), b 6 = 4.59509 (1 - i), 
a 7 = 0.118873 (1 - i), b 7 = 1.20969 (1 + i), a 8 = 0.118873 (1 + i), b s = 0.0511439 (1 - i) 

It is readily seen that all matrix elements of the pairing operator involving states 
'2' and '3' exhibit the square-root divergence if 5 — > 0, i.e. when approaching a 
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pseudo-DP. This divergence is cancelled out if one adds diagonal pairing matrix 
elements for these two states. However, off-diagonal matrix elements between 'reg- 
ular' states (states '1' and '4' in this case) and any combination of states '2' and '3' 
are divergent because: a 5 = a 6 * and 07 = a 8 *. 

Dependence of a real part of the pairing energy on the interaction strength g 
in the non-integrable pairing Hamiltonian (4) with g' / g = — 1/2 is plotted along 
the cut at (1Ze(g), -l/(4\/2)). Pairing energies in strongly mixed eigenstates '2' 
and '3' have opposite signs and diverge when a pseudo-DP is approached, i.e. if 
TZe(g) — > 0. On the other hand, the sum of pairing energies in states '2' and '3' 
is both positive and finite and changes smoothly close to a pseudo-DP. These fea- 
tures of pairing energy in the neighborhood of a pseudo-DP are essentially different 
from those discussed in connection with effects of an ordinary DP on pair-transfer 
amplitude [25] as the coupling between eigenvectors involved vanishes at the DP. 



6 Conclusions 

Besides accidental degeneracies, sharp level crossings in the hermitian eigenvalue 
problem are allowed only between states of different quantum numbers (differ- 
ent symmetries). For the non-hermitian (complex-symmetric) eigenvalue problem, 
such as found in OQSs, exact eigenvalue degeneracies appear also for states of the 
same quantum numbers/symmetries. These are: (i) the DP, when Riemann sheets on 
which live eigenvalues just touch each other [17], and (ii) the EP, when two Riemann 
sheets are entangled by the square-root type of singularity [20-22]. We have ana- 
lyzed here features of degeneracies in the non-hermitian and non-integrable 3-level 
pairing model (see Eq. (4) which is a schematic representation of the OQS Hamil- 
tonian. In particular, we have shown the appearance of a new kind of degeneracy, 
the pseudo-DP, in the anti-hermitian limit of this model as a result of the coales- 
cence between two EPs. The resonance eigenfunctions at the pseudo-DP (EP) are 
entangled and, except for the (complex) eigenenergy, it is impossible to define any 
physical quantity separately in each of the two resonances involved in the pseudo- 
DP degeneracy. 

This entanglement of wave functions may involve many states of an OQS, lead- 
ing to strong mixing of wave functions for states of largely different eigenener- 
gies. This mechanism is efficient not only in the close neighborhood of a particu- 
lar pseudo-DP (EP), i.e. for a particular choice of coupling constants in the OQS 
Hamiltonian, but also for systems with largely different coupling constants. Future 
studies, using the CSM and realistic effective interactions could provide a first at- 
tempt to provide a consistent description of continuum wave functions and related 
observables in the presence of a pseudo-DP (EP). 
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